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Abstract. Strangeness production in microscopic transport models for relativistic 
heavy-ion collisions from SIS to RHIC is reviewed: after a brief introduction into 
elementary strangeness production processes, the main emphasis is put on strangeness 
as indicator of the nuclear equation of state, the excitation function of the K + /ir + 
ratio and strangeness as a deconfinement indicator. 


1. Introduction and overview 

Microscopic transport models are a unique tool for the study of relativistic heavy-ion 
collisions: they offer a method of connecting the observable final state of such a collision 
with the time-evolution of the reaction, its dynamics and many not directly observable 
but very much sought after quantities and phenomena like the nuclear equation of 
state and the deconfinement phase-transition to a Quark-Gluon-Plasma [1]. Progress in 
heavy-ion physics requires a symbiosis between experiment and (transport-)theory. 

The basic principles of microscopic transport models can best be explained by 
studying the Vlasov-Uehling-Uhlenbeck (VUU) equation: 




x imi - /oil - h) - /i/ 2 (i - /;>(i - m 


This equation describes the time-evolution of the one-particle distribution function 
/i(ri,pi,t) which contains the positions and momenta of all particles of the system 
(i.e. the microscopic degrees of freedom) - initially all the protons and neutrons (or, in 


the case of a partonic description, the quarks and gluons) of the two colliding nuclei. 


These particles may interact either through an interaction given by the the real part of 
the retarded self-energy 9?E + (often also approximated by and referred to as the mean 
field) or through binary scattering, symbolized in this equation by the differential cross 
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Figure 1. Overview of available microscopic transport models and their range of 
applicability in incident beam energy. 


section If the r.h.s of this equation, the collision integral, is neglected one obtains 
the Vlasov equation. On the other hand, if the real part of the retarded self-energy !RS + 
is neglected and if the so called Pauli blocking factors (1 — /) in the collision integral 
are approximated by 1, one obtains the famous Boltzmann equation. 

In terms of the details of the implementation of mean held and collision term (e.g. in 
the number of hadronic resonances included or in the parametrization of cross sections), 
transport models may vary widely. Figure 1 provides an overview of currently available 
microscopic transport models applicable from SIS up to RHIC energies: models with 
purely hadronic degrees of freedom utilizing both, a mean held as well as a collision term 
are (I)QMD [2], VUU/BUU [3], ART [4] and BEM [5]. Models not containing a mean 
held are commonly referred to as cascade models - ARC [6] and LUZIFER [10] fall into 
this category. At higher energies initial particle production requires the introduction of 
string excitations - models which combine strings and hadrons are RQMD [7], UrQMD 
[8] and HSD [9]. At RHIC energies strings and hadrons may not anymore be the 
relevant degrees of freedom - here deconhnement needs to be taken into account and 
the elementary degrees of freedom for the initial reaction stage are quarks and gluons 
interacting through hard scattering as in VNI/BMS [11] and AMPT [12], 

2. Mechanisms of strangeness production 

Strangeness may be produced either in initial collisions among the incoming nucleons 
of the two colliding nuclei or through secondary interactions among produced particles, 
























Strangeness Production in Microscopic Transport Models 


3 




Figure 2. Left: excitation function of the exclusive p + p p + A + I\ + production 
cross section close to threshold in a resonance model. Right: K~ +p scattering cross 
section in UrQMD compared to data. 

e.g. pions and nucleons or excited resonances. 

2.1. Strangeness production at threshold 

There are two different approaches for initial strangeness production close to threshold: 

(i) direct parametrization : the individual cross sections for different exclusive 

production channels, e.g. p + p — > p + A + K + are directly parameterized. This 
method is very accurate if the cross sections are well known (i.e. well measured) but 
may become cumbersome away from threshold with increasing number of explicit 
channels to parameterize. 

(ii) resonance model here strangeness production is a two step process. Initially a 
heavy baryon resonance is excited, e.g. via p + p —>■ N + Nf 7l0 , which subsequently 
decays via A”* 710 —> Y + K + . This approach allows for an easier extension into the 
higher energy domain and may provide some rudimentary guidance for unknown 
strangeness production cross sections in secondary collisions, e.g. tt + N —> iV 1 * 710 —■» 
Y+K + . The left frame of figure 2 shows a fit of the resonance model to the exclusive 
p + p —> p + A + K + reaction channel. 

Secondary interactions like pion-induced strangeness production or flavor-exchange 
reactions are at least as important for the reaction dynamics and final strangeness yield 
as the initial/primordial strangeness production channels. The right frame of figure 2 
shows the K~ + N reaction cross section, which exhibits a distinct resonance structure. 
The hyperon resonances which are excited via this cross section may then either decay 
again into the K~ + N channel, or to almost equal probability decay into the Y + n 
channel, thus transferring strangeness in and out of baryonic degrees of freedom. As we 
shall later see, this K~ + N <-»•Y + tt exchange reaction is of particular importance. 
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Figure 3. Left: excitation function of the K and 7r multiplicity in elementary p + p 
reactions in UrQMD at large CM-energies compared to data. Right: distribution 
of string masses in Pb+Pb reactions at the CERN-SPS (160 GeV/u) and the fl~ 
production probability as a function of string-mass for two different values of the 
string-tension k. 


Ambiguities in both types of approaches arise from unknown and mostly 
immeasurable cross sections such as strangeness production in interactions involving 
baryon- and meson-resonances, e.g. A 1232 + N —* K + + X. 

2.2. Strangeness production at high energies 

At higher incident beam energies, particle production in general is dominated by string 
excitation and fragmentation. Up to 70% of the total strangeness produced in a Pb+Pb 
collision at top CERN-SPS energies is produced in initial highly energetic nucleon- 
nucleon interactions which lead to the excitation and subsequent fragmentation of strings 
[ 8 ]. On an elementary hadron-hadron level, the parameters of the string fragmentation 
are fitted to measured multiplicities and momentum distributions. The left frame of 
figure 3 shows the resulting excitation function of the pion and kaon multiplicity in 
proton-proton reactions in the UrQMD model (plot symbols are data, lines the model 
fit). 

When going from elementary hadron-hadron interactions to collisions of heavy 
nuclei one can explore the possibility of medium effects due to the hot and dense 
environment of the reaction. One such (non-hadronic) medium effect is the formation of 
color-ropes [13] - overlapping strings form a chromo-electric field which due to its larger 
field strength compared to the individual color-fluxtubes has an enhanced probability 
of fragmenting into strange hadrons. Rope effects can be simulated by increasing the 
string-tension k from it’s vacuum value of k = 1 GeV/fm to a value of k = 3 GeV/frn 
for strings fragmenting in the high energy-density region of the collision. The right 
frame of figure 3 shows the distribution of string masses excited in a Pb+Pb collision 
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Figure 4. Left: time-evolution of the K~ yield with and without inclusion of 
potentials. Right: ratio of energy-spectra of K + in Au+Au and C+C collisions for 
different cross section and equation of state parametrizations. 


at 160 GeV/u as well as the fragmentation probability of the string into an fl _ , both 
for k — 1 GeV/fm as well as for k = 3 GeV/fm [14]. A drastic enhancement of the fl - 
formation probability by more than one order of magnitude is visible. 

3. Strangeness and the nuclear equation of state 

In the SIS and AGS energy domain the main issues to be addressed by microscopic 
transport theory with respect to strangeness are the existence of in-medium 
modifications for kaons and the nuclear equation of state. In particular the K~ is 
considered to be strongly influenced by a scalar as well as a vector interaction. The left 
frame of figure 4 shows the time-evolution of the K~ yield for different combinations 
of potentials acting upon the K~ or the K + [15]. Surprisingly the K~ yield is totally 
insensitive to K~ potentials but very sensitive to K + potentials. This flavor coupling 
can be well understood through the flavor-exchange reaction K~ + N Y + n with the 
the hyperons being linked to the K + and their potential via associated production or 
higher resonance states like N* 7W Y + K + . The main observable to explore potential 
effects on K~ is the so called kaon-flow, which has been shown to be extremely sensitive 
to the strength and nature of the real part of the kaon-nucleon interaction [16]. 

The K + yield, on the other hand, has been found to be very sensitive to the 
nuclear equation of state [17]. Ambiguities may arise due to different implementations 
of unknown production cross sections like A 12 32 + N —> N + Y + K + . In order to 
avoid these, ratios of K + spectra for heavy and light collision systems can be used 

[18] . The right frame of figure 4 shows the ratio of the K + energy spectra for Au+Au 
and C+C collisions for different cross section and equation of state parametrizations 

[19] . A strong sensitivity to the equation of state is observed while different production 
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Figure 5. Excitation function of the K + /n + ratio in RQMD (left) and HSD (right; 
open symbols refer to HSD. closed symbols are data). 


cross sections do not alter the value or functional form of the ratio. Comparisons of 
current calculations with data favor a soft equation of state with momentum-dependent 
interaction [18]. 

4. Excitation functions 

The K + /tt + ratio provides a measure of the ratio of newly produced strange to non- 
strange valence-quarks in a heavy-ion reaction. A kink in its excitation function vs. 
incident beam energy was thought to hint at a possible deconfinement phase-transition. 
Recent calculations in the framework of a statistical model show that the K + /n + ratio 
is expected to reach its maximum value around a beam energy of 30 GeV/u, well in 
line with experimental hirelings [20]. While the agreement of the data with statistical 
model calculations cannot per se prove or disprove the existence of a deconknenrent 
phase-transition, the K + /tt + ratio as deconfinement indicator remains ambiguous at 
best. 

However, the beam-energy dependence of the K + /n + ratio has emerged as a rather 
stringent test for microscopic transport models: figures 5 and 6 show the excitation 
function of the K + /7t + ratio for RQMD (fig. 5 left, provided by [22]), HSD (fig 5 right, 
taken from [23]) and UrQMD (fig. 6 left, provided by [21, 22]). While RQMD compares 
well to the data, both HSD and UrQMD exhibit problems which may not necessarily 
be rooted explicitly in the strangeness production mechanisms of the particular models: 
in the case of HSD, an analysis shows that the observed functional dependence of the 
K + / 7T + ratio is very sensitive to the threshold of the string production cross section and 
the excitation of the high-mass resonance continuum [24], Older versions of RQMD with 
a similar implementation of these cross sections exhibit the same behavior as observed 
in HSD. In the case of UrQMD the under-prediction of the K + /tt + ratio starting at AGS 
energies can be traced to a pion excess (the kaon yields roughly agree with the data). 
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Figure 6. Left: excitation function of the K + /tt + in UrQMD. Right: comparison of 
the pion per participant excess in Pb+Pb reactions vs. p+p reactions in UrQMD to 
data. 

The right frame of figure 6 shows the pion per participant excess in Pb+Pb reactions 
vs. p+p reactions in UrQMD compared to data [21, 22]: in UrQMD this excess sets in 
too early with respect to the incident beam energy. It is very likely that this behavior 
is caused by the violation of detailed balance due to the lack of multi-particle (mostly 
pion) collisions which would convert the pion excess into heavier particles such as anti- 
baryons and hyperons [25, 26] - preliminary calculations incorporating these processes 
yield a Ii + /ir + ratio very close to the data [21]. 

A more detailed study of various excitation functions of the K + /n + ratio in RQMD 
as well as a discussion of the K + /K~ ratio can be found in [27]. 

5. Strangeness as deconfinement indicator 

The relative enhancement of strange and especially multistrange baryons with respect 
to peripheral (or proton induced) interactions has been suggested as a signature for 
the transient existence of a QGP-phase [28, 29, 30]: the main argument being that the 
(chemical or flavor) equilibration times should be much shorter in the plasma phase 
than in a thermally equilibrated hadronic fireball of T ~ 160 MeV. 

The dominant production mechanism in an equilibrated (gluon rich) plasma phase, 
namely the production of ss pairs via gluon fusion (gg —> ss) [28], should allow for 
equilibration times similar to the interaction time of the colliding nuclei, and to the 
expected plasma lifetime (a few fm/c). 

The yields of strange baryons per event calculated in UrQMD are shown in the 
left frame of figure 7 as a function of the number of participants for Pb+Pb and 
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Figure 7. Left: excitation function of strange baryon multiplicity vs. number of 
participants in UrQMD. Right: Inverse slopes of the tut -spectra of 7r, K. p, A + S°, 
S° + S”, and Cl~ at y c . m . = 0. 

p+Pb collisions at 160 GeV/u [31]. The A + A- (circles), 5“ + S - - (squares), and 
0 — + 0 — - (triangles) values are shown. The stars correspond to experimental data of the 
WA97 collaboration [32], Open symbols represent the results of the standard UrQMD 
calculations, whereas full symbols exhibit a calculation with an enhanced string tension 
of k — 3GeV/fm, for the most central collisions (A^ part > 300). Obviously the standard 
UrQMD calculation, which can be seen as a baseline of known hadronic physics, strongly 
underestimates the (multi-)strange particle yields in central collisions, in particular in 
the case of the fU. Only the inclusion of non-hadronic medium effects, like color-ropes 
[13], which are simulated by increasing the string-tension for central collisions (see also 
the right frame of figure 3), enhances the yield to a level of near-compatibility with the 
data. Similar hirelings have also been made in the context of HIJING calculations [33]. 
While these findings by no means prove the validity of the color-rope approach, they 
clearly show the necessity of some kind medium effect beyond regular binary hadronic 
(re)scattering in order to understand the data. This statement is corroborated by a 
calculation with the string model NeXuS, including so-called QGP droplets - domains 
of high energy-density hadronizing according to a statistical phase-space population - 
which give rise to a similar strangeness enhancement [34], 

However, recently hadronic multi-particle interactions in the early dense reaction 
phase have been suggested to significantly enhance the yield of anti-protons and (anti- 
)hyperons [25, 26]. It remains to be seen, however, whether these effects are sufficient to 
explain the observed enhancement or other non-hadronic (i.e. deconfinement based) 
effects need to be taken into account. 

The study of deconfinement and a subsequent phase-transition to deconhned 
hadronic matter poses a great challenge to microscopic transport models. In most 
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approaches hadronization is a uni-directional process and the equation of state of the 
system ill-defined. One possible remedy is to use a hydrodynamical approach for 
the early deconhned phase of the reaction and subsequent phase-transition, coupled 
with a microscopic calculation for the later, hadronic, reaction phase in which the 
hydrodynamical assumptions are not valid any longer [35]. In the following, such a 
combined macro+micro model will be used to study the flavor- and mass-dependence 
of hadronic slope parameters: the right frame of fig. 7 displays these inverse slope 
parameters T obtained by an exponential fit to dNi/d 2 rriTdy in the range rrir — rii l < 
1 GeV for SPS and RHIC [36] and compares them to SPS data [32], The trend of the 
data, namely the “softer” spectra of S’s and Q’s as compared to a linear T(m ) relation 
is reproduced reasonably well. This is in contrast to “pure” hydrodynamics with kinetic 
freeze-out on a common hypersurface (e.g. the T = 130 MeV isotherm), where the 
stiffness of the spectra increases linearly with mass as denoted by the lines in fig. 7. 
When going from SPS to RHIC energy, such a hybrid model as discussed here generally 
yields only a slight increase of the inverse slopes, although the specific entropy is larger 
by a factor of 4-5 ! The reason for this behavior is the first-order phase transition that 
softens the transverse expansion considerably. 

The reason for the softening of the spectra is that the hadron gas emerging from 
the hadronization of the QGP is almost “transparent” for the multiple strange baryons. 
Analyzing the collision numbers, one finds that Q’s suffer on average only one hadronic 
interaction, whereas TV’s and A’s suffer far more collisions with other hadrons before 
they freeze-out. Thus, one may conclude that the spectra of S’s and especially Q’s 
are practically unaffected by the hadronic reaction stage and closely resemble those 
on the phase boundary. They therefore act as probes of the QGP expansion prior to 
hadronization and can be used to measure the expansion rate of the deconhned phase. 
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